Adsorption Transition of a Polyelectrolyte on a High-dielectric Charged Substrate 
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The behavior of a polyelectrolyte adsorbed on a charged surface of high-dielectric constant is stud- 
ied by both Monte-Carlo simulation and analytical methods. It is found that in a low ionic strength 
medium, the transition is first-order with the repulsive charged surface. The surface monomer den- 
sity, which is the order parameter of the adsorption transition, follows a linear relation with surface 
charge density. It indicates the polyelectrolyte is compressed on the substrate without any confor- 
mational change before the de-sorption. Finally, a new scaling law for the layer thickness is derived 
and verified by simulation. 

PACS numbers: 61.25.Hq, 82.35.Gh 



Polymer adsorption on an attractive surface has drawn 
considerable interest due to its relation to surface effects 
in critical phenomena and practical importance in ma- 
terial science and biophysics. It is well established that 
the adsorption transition is continuous if its attraction 
on the surface is short-ranged 0, 0> 0, 0] • On the other- 
hand, long-ranged electrostatic interactions in polyelec- 
trolyte systems pose many challenging theoretical prob- 
lems. Recently the macroion adsorption on an electro- 
statically attractive interface and the associated charge 
inversion phenomena of adsorbed polyelectrolytes acquire 
lots of attention 0,0,0- 

Previous analytical approaches using the Edwards 
equation imposed the continuity of the monomer den- 
sity across the surface and setting the monomer density 
to zero at the surface 0, Q . Within the framework of the 
self-consistent field method, both the Poisson-Boltzmann 
equation and Edwards equation were solved simultane- 
ously 0, 0, 0| with zero monomer density at the sur- 
face. These treatment however cannot faithfully respect 
the electrostatic boundary condition. As a result, the 
adsorption transition would still be continuous whatever 
the surface potential looks like. 



Recently, the theoretical interest to the problem is due 
to its importance for multi-layer polyelectrolyte adsorp- 
tion yjj LU3| • It also raises the question of applying 
Poisson-Boltzmann theory to polyelectrolyte adsorption 
because the theory fails to capture the correlation effects. 

In this letter, we study the adsorption of a single poly- 
electrolyte on a high-dielectric substrate in which the im- 
age charge attraction is strong. At low ionic strength, the 
adsorption transition occurs when the surface charges are 
repulsive instead of the attractive case that were usually 
studied. The problem is tackled by performing Monte- 
Carlo simulations and also by analytical methods in poly- 
mer physics taking full account of the appropriate bound- 
ary conditions. It is found that the order of the adsorp- 
tion transition, the physical mechanism, and the scaling 
behavior are all different from those of the attractive sur- 
faces. 

A polyelectrolyte carrying positive charges is immersed 
in a medium [z > 0) of dielectric constant e. At z = 
there is an impenetrable surface of uniform surface charge 
density a. Below the surface (z < 0), it is a substrate of 
dielectric constant e'. Denote the charge on a polymer 
segment ds by qods, the Hamiltonian is 
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where s is the variable to 



parametrize 
l 



the 



chain, 

Iq the bare persistence length, and the Debye 

screening length. f(s) = (x(s),y(s), z(s)), f'(s') = 
(x(s'),y(s'), —z(s')) are the positions of the monomers 
and their electrostatic images, respectively. T = <?g/e, 
T' = r(e' - e)/(e' + e), and h = 47rq er/(e' + e) are the 



coupling parameters governing the strength of Coulomb 
interactions among the monomers themselves, between 
the polymer and its image, and between the polymer and 
the charged surface, respectively. The last term in Eq.QJ 
represents the excluded volume interactions with u> > 
(good solvent regime) in this study. We shall focus on the 
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case of a charge polymer in a low ionic strength medium. 

The above continuum model is discretized to perform 
Monte-Carlo simulation. The continuous curve f(s) is 
replaced by a chain of beads r*j (i — 1, . . . , N) with hard- 
core excluded volume of finite radius a. Total lengths 
up to N — 120 are employed. Units of length and en- 
ergy are set to be 2a and q^jlta, respectively. Dielectric 
ratios e'/e are chosen from 2 to 12.5 (aqueous solution 
with a metallic substrate) . Runs up to 10 9 MC steps are 
performed to achieve good statistics. 

The adsorption layer can be characterized by the nor- 
malized monomer density p(z). p a = p(a) representing 
the fraction of monomers being adsorbed on the substrate 
is chosen as an order parameter to describe the adsorp- 
tion transition. p a > and p a = characterize the 
adsorbed and de-sorbed states respectively. In FigQ] p a 
as a function of the surface charge density a for various 
e'/e > 1 is shown. The discontinuous jump of p a across 
the threshold indicates the transition is first order. We 
also verified that the energy jump (latent heat) across 
the transition is proportional to N. Similar results were 
obtained for larger 
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FIG. 1: Monte-Carlo results for the normalized monomer den- 
sity at the surface, p , as a function of surface charge density 
a (in unit of qo/Aa 2 ) for different e'/e at k" 1 = 25. The fit- 
ted straight lines are terminated at their adsorption transition 
points. The vertical dashed lines are drawn as guides to the 
eyes. Inset: The polarization surface charge density induced 
by the polyelectrolyte, |<r p |, as a function of dielectric ratio 
e'/e. The sign of a p is opposite to qo and is negative. The 
solid curve is fitted from Eq.JHJ with cr po i y = 0.118. 

Furthermore, the data in Fig^ also indicate that p a 
is linear in a with the slope depending on the ratio of 
e'/e. Such a linear relation between p a and a can be 
understood from the electrostatic boundary conditions 
that the system has to satisfy. The electric potential 
4>{z) in the neighborhood of the z = boundary obeys, 
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where er p is the polarization surface charge density in- 
duced by the polymer only, which depends on e'/e but is 
independent of a in the adsorbed regime near the tran- 
sition. Notice that er p in general is a complicated func- 
tion since it relies on the polymer conformation. Also, if 
one treats the polymer as a marcomolecule with a well- 
defined surface, its surface charge density at z — a should 
be proportional to the monomer density p a , this also ap- 
plies to the electric field in the z < region, 
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where K > is the corresponding proportional constant. 
Applying the continuity condition for the electric field 
from z = + to z = a~ , and using Eqs.(0|-l@J), one gets 
the linear behavior 
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Notice that K and a p are functions of e'/e. The linear 
behavior in Eq. (JSJ is confirmed by the simulation data 
as shown in Fig^ the slope decreases monotonically with 
e'/e. Substituting a = into Eq.JSJ), we get the polar- 
ization surface charge density as a function of dielectric 
constant ratio, 
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K for different dielectric ratios is obtained from the slopes 
of different straight lines presented in Fig.l. a p as a 
function of e'/e is then fully determined from simulation 
data as shown in the inset of Fig.l. It fits very well 
to Eq.JfjJ with cr po i y = p a \o=o/K = 0.118. It suggests 
the conformation of the adsorbed polymer is compressed 
on the high-dielectric substrate. We have also check the 
adsorption layer thickness is independent of the number 
of monomer N, which is consistent with the picture of 
a compressed state. There is no conformational change 
during the de-sorption contrary to the scaling results as 
predicted by Borisov et. al. 0] for the attractive surface. 

The polyelectrolyte behaves as electric blobs arranged 
longitudinally and lie down parallel to the surface. In- 
creasing the attraction from the charged surface reduces 
the ^-fluctuation amplitude of the chain, but the effective 
in-plane surface charge distribution of the polyelectrolyte 
does not change. Hence the polarization a p is indepen- 
dent of a. The excluded volume effect is safely ignored 
because it takes almost no effect in the z-direction. The 
effect from self-electrostatic interaction of the polyelec- 
trolyte can be absorbed into the bare persistence length 
from Iq to I. 

Because the monomer would feel the strongest attrac- 
tion from its direct image around the adsorption regime, 
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the r'-term in Eq. (JTJ is approximated by the interaction 
of each monomer and its corresponding image only. The 
residual attraction from the images of other monomers is 
then adsorbed into the coupling parameter V from q to 
q. The partition function is reduced to 
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Transforming the variable from r(s) to the normalized 
monomer density p(r) = dsS(f — f(s)) by in- 

troducing an auxiliary field, then applying the ground 
state dominance approximation in large- iV limit and by 
variational principle 0, 0] , one obtains the Edwards- 
Schrodinger equation, 
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where £o acts as a Lagrange multiplier to enforce the con- 
straint of the ground state wavefunction normalization. 
The monomer density is given by p(z) = \ip(z)\ 2 . Eq.JSJ) 
also describes a quantum particle at its ground state 
moving under a combined potential of a 1-d screened 
Coulomb attraction and an almost linear potential. How- 
ever, the boundary condition expressed by Eq.© is dif- 
ferent from the hard- wall boundary condition ip\ s = 
usually employed for a quantum particle. Instead ip\s = 
^fp~a~ 7^ for the present problem implies that the steric 
force felt by the polyelectrolyte from the charged surface 
should be modified Setting ip\ a = HIIIIII in the 
polyelectrolyte adsorption problems in previous studies 
are not completely correct. 

During the adsorption, the rod-like polyelectrolyte 
tends to lie down on the charged surface. The thick- 
ness of the adsorption layer is of the same order of the 
gyration radius in z-direction. At low ionic strength in 
which the Debye length is much greater than the layer 
thickness, the polyelectrolyte cannot feel the potential of 
length scale much larger than but only the potential 
barrier height is important. The original potential V(z) 
in Eq.lJBJ can thus be replaced by 
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where Zbr is chosen such that V'(zb r ) — and V(zb T ) is 
the barrier height. In the limit of cr = K = 0, analytic 
solution gives 
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where W>, : i/2 is the Whittaker's notation of the confluent 
hypergeometric function [2^, and A is the least value 



satisfying the boundary condition. Bound state exists 
for arbitrary e'/e > 1. It implies the threshold surface 
charge density, a t > at low ionic strength. 

For both a, n > 0, no exact solution exists in general 
but one can analyze it around the transition. Near the 
surface, the image charge attraction dominates over the 
surface charge repulsion and hence the binding energy is 
approximated by the a — case in Ea. lfTtHl . The poly- 
electrolyte undergoes a de-sorption transition when the 
binding energy meets the barrier height V^Zbr). After 
some algebra, we have at ~ (e'/ e — 1) f° r e '/ e ^ 1 an( i 
ft ~ ( e '/ e — I) 3 f° r e 'A ?i 1- This analytic result is 
consistent with our simulation data as shown in Fie[2Ji. 
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FIG. 2: (a) Simulation results for the surface charge density 
at the transition, at, as a function of dielectric ratio e'/e in 
logarithmic scale at k _1 = 25. The straight lines indicate 
slopes of 1 and 3 as suggested in the text, (b) Simulation 
results for the inverse decay length a (which is proportional 
to the inverse layer thickness) as a function of a (in unit of 
qo/4o 2 ) for e'/e = 12.5, ft -1 = 25. The straight line is a linear 
fit. <Tt = 0.102 in this case, a is obtained from exponential 
fitting to the tail of corresponding density profile. 

An approximate solution for the density profile p{z) 
for the a > case can be obtained by variational method 
with trial wavefunction 



i>{z) = + fia(z - a))e-^ z - a > 
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where a -1 is the decay length, /i is positive because the 
trial wavefunction is restricted to be nodeless. a and 
/i are not independent but related via the wavefunction 
normalization condition. The inverse decay length is cal- 
culated to be 



a = 3/3T'/2l 2 
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where the leading term is independent of a. Near the 
transition, the decay length and hence the thickness of 
the adsorption layer increase and remain finite. From 
Eq.lJSJ and (|12|) . we get the scaling behavior 



a — at ~ (at — a) for < a < a t 



(13) 
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where at is the threshold inverse layer thickness at a = 
(7,7 . The variation of a as a function of a obtained from 
the simulations is shown in Fig|2) which can be well fit- 
ted to a linear relation consistent with Ea. (|13l) . On the 
other hand, for the case of adsorption onto an attractive 
charged surface (<r < 0) with substrate of e'/e < 1 (e.g. 
DNA in aqueous solution adsorbed onto a charged lipid 
membrane) , the asymptotic solution to Eq. JHJ) reproduces 
the usual scaling a ~ |cr| a and is a continuous adsorption 
transition 0] , and the thickness swells to infinity as the 
polyelectrolyte is de-sorbed. 

A strongly charged polyelectrolyte immersed in a salt 
solution will attract oppositely charged ions to condense 
until its effective charge density reaches the Manning 
threshold [13 • This means that one can just renormalize 
<7o in our system to 2ea/l^ if qo is larger than 2ea/l^ (7b 
is the Bjerrum length). Similarly, the strongly charged 
surface of bare charge density larg er than k/ (ttIb) is just 
renormalized back to k/(ttI-b) |24j . The Gouy-Chapman 
length, calculated from the nonlinear Poisson-Boltzmann 
theory, is of the order (Ibt)^ 1 , which is very large 
around the transition. However, the nonlinear Poisson- 
Boltzmann potential [25j near the substrate is given by 
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where Aqc is the Gouy-Chapman length and 7 = 
y/ k 2 Aq C + 1 — ft Ago Notice that the linear term in z 
is proportional to AqJ, oc a and is identical to the lin- 
ear term as expanded from the Debye-Hiickel potential. 
It is not surprising since both Poisson-Boltzmann and 
Debye-Hiickel potentials share the same boundary condi- 
tion. When the adsorbed polyelectrolyte layer thickness 
is of one to two monomer size (as seen from a -1 ~ 1 in 
Fig. 2b), the surface potential felt by the polyelectrolyte 
should be linear. Physically speaking, it does not matter 
whether the potential is Poisson-Boltzmann or Debye- 
Hiickel, or even the linear one if the layer thickness is 
much smaller than Agc and n^ 1 . The potential near 
the surface is determined from the boundary condition. 
Unlike the case of an attractive surface |26(, the effect 
from large Guoy-Chapman length near the transition is 
irrelevant in the high-dielectric case. 

Our results on the single polyelectrolyte adsorption 
may provide a starting point to study the charge in- 
version and multi-layer adsorption 27]. At low ionic 
strength, polyelectrolytes are adsorbed in a multi-layer 
structure because of strong Coulomb repulsion. Each 
layer is composed of parallel Id Wigner crystal [28j |. 
The upper bound of the multi-layer thickness is zt, r ~ 
ct~2 (e'/ e — 1)2. It suggests we can easily adjust a single 
layer adsorbed onto a high-dielectric substrate by tuning 
the surface charge density. Rigorous treatment based on 
this physical picture will be elaborated elsewhere. 

In conclusion, the adsorption transition of a single 
polyelectrolyte on a high-dielectric substrate is first or- 



der since the polyelectrolyte needs to overcome a bind- 
ing energy from its image charge. Because of the strong 
Coulomb attraction as compared to the linear repulsive 
potential near the surface, the polyelectrolyte is com- 
pressed without any conformational change before the 
de-sorption. A scaling law for the adsorption layer thick- 
ness is also derived and verified by simulation. 
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